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ABSTRACT 

PSR J2032+4127 is a 7 -ray and radio-emitting pulsar which has been regarded 
as a young luminous isolated neutron star. However, its recent spin-down rate 
has extraordinarily increased by a factor of two. We present evidence that this 
is due to its motion as a member of a highly-eccentric binary system with a 
~ 15-Mq Be star, MT91 213. Timing observations show that, not only are the 
positions of the two stars coincident within 0.4", but timing models of binary 
motion of the pulsar fit the data much better than a model of a young isolated 
pulsar. MT91 213, and hence the pulsar, lie in the Cyg OB2 stellar association, 
which is at a distance of only 1.4-1.7 kpc. The pulsar is currently on the near 
side of, and accelerating towards, the Be star, with an orbital period of 20-30 
years. The next periastron is well-constrained to occur in early 2018, providing 
an opportunity to observe enhanced high-energy emission as seen in other Be-star 
binary systems. 
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1 INTRODUCTION 

PSR J2032+4127 is a 143-ms pulsar discovered by the 
Large Area Tele scope (L AT) of t he Fermi Gamma-ray 
Space Telescope (lAbdo et al.l 2009ll and subseque ntly de¬ 
tected at radio frequencies I Camilo et al.l FoOOl) . Using 
7-ray data as well as radio data, pulse arrival-time ana l- 
ysis by those authors, and latterly bv lRav et al.l d201 il l . 
showed that the pulsar had a large slow-down rate, in¬ 
dicating that it was a young pulsar, a notion which has 
been supported by a recent glitch in its rotation rate. 
The timing analysis showed that the projected position 
of the pulsar lay close to a V=11.95 Be star, MT91 213 in 
the C yg OB2 stellar association. However, ICamilo et al.l 
(l2009l l concluded that it was not a binary companion of 
PSR J2032+4127, based mainly upon the apparent lack 
of variations in the pulsar rotation rate that would arise 
from the Doppler effects of any reasonable (circular) or¬ 
bital binary motion. Subsequent studies of the proper¬ 
ties of the pulsar have all been based upon the assump¬ 
tion that it is a solitary young energetic pulsar. For in¬ 
stance, lAliu et ahl (l2014l l present a VERITAS detection 
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of TeV J2032+4130 and they discuss its likely associ¬ 
ation with PSR J2032+4127, noting that the extended 
nature of the source and the prevalence of pulsar wind 
nebulae (PWNe) in the Galactic TeV source population 
argues that it is a PWN powered by the pulsar. However 
they note that the extended X-ray e mission that is also 
spatially coincident with the source Jllorns et ah 1 120071 : 
iButt et al.l 120061 : iMurakami et al.ll201lh is quite weak if 
it is from a PWN. Six years of both 7-ray and radio tim¬ 
ing data are now available, and we revisit the possibility 
that the pulsar is in orbit with MT91 213. 


2 OBSERVATIONS AND TIMING NOISE 
MODEL 

Fermi LAT data were used from shortly after the 
start of the mission in August 2008 up to June 2014 
(MJD 5468 2-56824). F ollowing the maximum likelihood 
method of I Ray et 5] (120111 1 , we constructed times-of- 
arrival (TOAs) with a 14-day cadence from “reprocessed” 
Pass 7 Fermi LAT data. For this analysis, we used 
SOURCE class photons, excluding events with a zenith an¬ 
gle > 100° and when the spacecraft rocking angle exceeds 
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52°. To improve sensitivity, we employed photon weight¬ 
ing using the spectral models available in The Second 
Fer mi Large Area Tele scope Catalog of Gamma-ray Pul¬ 
sars (lAbdo et al.ll2013lh 

Radio timing observations were made with the 
NRAO Green Bank Telescope (GBT) and the Lovell Tele¬ 
scope (LT) at Jodrell Bank. The GBT observations were 
primarily in bands centred on 820 MHz and 2000 MHz 
during the early pa rt of this period (MJD 54836-55589) 
dCamilo et al .1120091) . as well as MJD 56855-56857. Radio 
observations with the LT were made at approximately 
weekly intervals at around 1520 MHz for MJD 55222- 
56830. As described bv ICamilo et al.l (I2009I ). the 7 -ray 
and radio pulse profiles are very different in form. In order 
to establish the true alignment of the profiles, the disper¬ 
sion measure (DM) was initially determined using only 
radio data. Inspection of the radio pulse profiles shows 
that there is little change in shape between 800 MHz 
and 2000 MHz, and a single standard profile was used 
to obtain TOAs for all the radio observations. The DM 
was determined at three epochs, using the three 800-MHz 
TOAs and their neighbouring 1520-MHz and 2000-MHz 
TOAs. This allowed the time offset of the contemporane¬ 
ous 7 -ray T OAs to be establish ed using the jump facility 
in tempo2 ((Hobbs et al.l l 200^ 1. The mean value of these 
ti me offsets was consis tent with the alignment obtained 
bv ICamilo et al.l (120091 ) and was applied to all the 7 -ray 
TOAs, allowing the delay between these TOAs and the 
radio TOAs to be used to determine the DM throughout 
the data set. 

The radio TOAs have errors of ~ 300 fis, about half 
the errors of those provided by the Fermi LAT. I 11 total, 
about 400 good TOAs are available over six years (MJD 
54682-56857). 

The pulsar shows substantial variations in rotation 
rate, indicating a large amount of “timing noise”, which 
is not unexpected in a young pulsar (i.e. having a large in¬ 
trinsic spin-down rate). However, it has recently become 
clear that the magnitude of these variations far exceeds 
any previous ly enc ountered timing noise in a pulsar (e.g. 
iHobbs et alj|20ld ). Fig. la presents the evolution of the 
barycentric rotation frequency of the pulsar over the 6 - 
year span of the timing data, showing a substantial re¬ 
duction in frequency during this time. Also apparent is 
an increase in the magnitude of the slope, corresponding 
to an increase in the magnitude of the frequency deriva¬ 
tive, shown in Fig. lb, by more than a factor of two, and 
a decrease in the characteristic age by the same factor, 
confounding any interpretation in terms of conventional 
pulsar spin-down. There is also a barely-discernable mi¬ 
nor glitch which occurred in Sept 2011 (MJD 55811) and 
which can be fitted by steps in frequency (A u) and fre¬ 
quency derivative (Ai>). The whole set of TOAs can be 
well-fitted by a model involving the position, three glitch 
parameters, the rotation frequency and 7 derivatives (Ta¬ 
ble 1, column 2). The timing residuals (the difference be¬ 
tween the observed and model-based TOAs) relative to 
this model, which we will henceforth refer to as the “noise 
model”, are shown in Fig. lc. 

The behaviour in the rotation seen in Figs, la and lb 
has never been seen in any other isolated pulsar, and is 
not characteristic of either a post-glitch recovery, during 
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Figure 1. The spin-frequency history of PSR J2032+4127 
determined over six years, a) The observed variation in the 
spin-frequency offset A v from 6.98083 Hz, determined from 
fits to 150-day sets of TOAs every 50 days, showing the mono¬ 
tonic spin-down of the pulsar. This spin-down is interrupted 
briefly by a small glitch close to MJD 55811, identified by 
an arrow, b) The variation of the spin-frequency first deriva¬ 
tive i> determined over the same time intervals as (a), show¬ 
ing a doubling of the magnitude. The fits for v and 0 for 
(a) and (b) were conducted simultaneously with the position 
fixed at a = 20 h 32 m 13 s .105, <5 = 41°27'24",36. c) The timing 
residuals relative to the 7-derivative noise model presented in 
Table 1, column 2. d) The timing residuals relative to a 6- 
derivative noise model, e) The timing residuals relative to the 
best-fit “binary” model, which has orbital period / : 7=8578 d 
and mass function / m = 10 Mg (Table 1, column 4). In pan¬ 
els c, d and e, the colour of the points indicates whether they 
were obtained from the LAT (red), the LT (black), the GBT 
at 2000 MHz (green) or the GBT at ~ 820 MHz (blue). All 3 
timing models used in Figs lc, Id, and le do in fact include 
fits for the glitch parameters. In the timing-noise models, there 
is significant covariance between glitch parameters and other 
polynomial parameters, so that some features of the glitch are 
absorbed in the polynomial parameters, leaving some sharp, 
unmodeled features of the glitch in the residuals. It seems that 
the glitch and Keplerian parameters are not so covariant, and 
the glitch is fitted well. 
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Figure 2. The orbital variation of the frequency v and the first 
derivative i> of PSR J2032+4127 for the best-fit binary model, 
which has orbital period P ^ = 8578 d and mass function / m = 
10 Mq (Table 1, column 4). The scale on the right-hand axis of 
the upper panel indicates the magnitude of the orbital radial 
velocity. The dotted line in the lower panel shows the variation 
in v on a scale expanded by a factor of 100. The bold sections 
in the two plots indicate the range of the available data. 


which the magnitude o f the slow-down rate is usually ob¬ 
served to decrease (e.g. [Espinoza et al.l201lh , or of timing 
noise, which is usually characterised by switches, often 
quasi-periodic, between values of slow -down rate rather 
than the s mooth variation seen here (jLvne et al.1 l2010l : 
lLvneil2013h . Most anomalous of all is the doubling of the 
slow-down rate, also not seen in any other isolated pulsar. 


3 A BINARY MODEL 

We believe that the only plausible origin of such a large 
variation in the observed slow-down rate of a pulsar lies 
in the Doppler effects of binary motion wit h another star. 
While we note the remark by ICamilo et all J2009h that 
there is no evidence of any short-period binary motion 
and that the binary period (of any circular orbit) must 
be in excess of 100 years, we have explored the possibility 


that the pulsar is actually a member of a long-period 
binary system with a large orbital eccentricity. 

First, we sought fits of binary models to the rotation- 
frequency data of Fig. la. While good fits to the data were 
possible using models for eccentric binary orbits with or¬ 
bital periods Pb in excess of about 6000 days (16 years), 
it soon became clear that there were strong covariances 
between some of the fitted parameters, arising from the 
small orbital phase range of the available data. In partic¬ 
ular, there were large covariances between the intrinsic 
pulsar slow-down rate, v, the orbital period Pb and the 
projected semi-major axis of the orbit x = (a/c) sin i, 
where a is the semi-major axis of the orbit, i is the in¬ 
clination of the plane of the orbit to the plane of the 
sky and c is the speed of light. These covariances allowed 
many good, but not unique, fits to the data. We therefore 
explored fits to the data for a range of fixed values of Pb 
between 6000 and 12000 days at 500-day intervals, and, 
for each of these, a range of fixed values of x, correspond¬ 
ing to fixed values of mass function / m of 2, 5, 10 and 20 
Mg. / m is a function of the masses of the neutron star 
(M p ) and its companion (M c ) and orbital inclination i 
and is determined from Pb and x from Kepler’s laws by: 

_ (M c sin i) 3 _ An 2 x 3 ,, 

- (M p + M c )2 ~ ~GPf 1 ’ 

where G is Newton’s gravitational constant. If M p , M c 
and / m are in solar masses, x is in light-sec and Pb in 
days, then 

x = 9.766 y fmP 2 ■ (2) 

Fig. 2 shows an example of one of the best fits to the 
data, for Pb = 8578 days, / m = 10 Mg. Note that the 
data span occupies only about 20% of the orbit. Remark¬ 
ably, the root mean square (rms) of the frequency resid¬ 
uals (the differences between the measured and model 
frequency values) for this simple model was only approx¬ 
imately 1 CP 5 of the total frequency variation during this 
time and is consistent with the measurement errors. 

For each Pb, /m pair, the ephemeris resulting from 
the frequency fit was subsequently used as the basis of 
a coherent timing analysis of the TO As using tempo2. 
Keeping Pb and x at their fixed values, the three other 
Keplerian parameters, the pulsar position, DM and its 
first derivative, DM1, the rotation frequency and its first 
derivative and three glitch parameters were all fitted to 
the TOAs. Fig. 3 summarises the results of the fits. In 
particular, Fig. 3a shows the rms of the timing residuals 
relative to the fitted model as a function of Pb and / m . 
There is a broad minimum in the rms with a value of 
about 0.6 ms for 7500 d < Pb < 9500 d and all values of 
/ m > 5 Mg provide indistinguishably good fits. Best-fit 
models were obtained for / m = 10 Mg and / m = 20 Mg 
by including a fit for Pb and the results are given in Ta¬ 
ble 1, columns 4 and 5. Because of the strong covariance 
between some of the parameters, we used a Monte-Carlo 
method to estimate the errors in the fitted parameters 
that are given in the table. The Monte-Carlo analysis 
used simulations of the pulsar observations seeded with 
the best-fit parameter file. The simulated data sets had 
the same epochs of the combined LAT and radio TOAs, 
giving the same cadence and orbital coverage as the real 
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Table 1 . Model fits to the TOAs of PSR J2032+4127. Observed values of parameters for a 7-derivative polynomial “Noise Model” 
fit, and for two best-fit “Binary Models”, constrained to have mass functions / m = 10 and 20 Mq respectively and with errors 
determined from a Monte-Carlo analysis. One-cr uncertainties are given in parentheses after the values and are in units of the least 
significant quoted digit. 


Parameter 

Noise Model 

Parameter 

Binary / m = lOM© 

Binary / m = 20M© 

Right Ascension, a (J2000.0) 

20 h 32 m 13. s lll(8) 


20 h 32 m 13. s 107 (8) 

20 h 32 m 13. s 105(8) 

Declination, S (J2000.0) 

41°27'24".17(10) 


41°27'24".36(10) 

41°27'24".36(10) 

Epoch of frequency, to (MJD) 

55714.0 


55700.0 

55700.0 

Freq, vq (Hz) 

6.9808479486(10) 


6.980808(12) 

6.980806(12) 

Freq 1st deriv, z>o (10 -12 s -2 ) 

-1.16583(12) 


-0.61(5) 

-0.61(5) 

Glitch Epoch, T g (MJD) 

55810.51(2) 


55810.76(7) 

55810.76(7) 

Freq, A u g (10~ 6 Hz) 

1.907(2) 


1.907(2) 

1.907(2) 

Freq 1st deriv, Ai> g (10 -16 s -2 ) 

-4(3) 


-10(2) 

-10(2) 

DM (pc cm -3 ) 

114.66(3) 


114.65(3) 

114.65(3) 

DM deriv, DM1 (pc cm -3 y -1 ) 

Timing noise parameters: 

0.00(1) 

Binary parameters: 

0.00(1) 

0.00(1) 

Freq 2nd deriv, o (10 -20 s -3 ) 

-0.3651(7) 

Orbital period, P\> (d) 

8578(1200) 

8625(1200) 

Freq 3rd deriv, (10 -28 s -4 ) 

-0.439(2) 

Epoch of Periastron, To (MJD) 

58161(31) 

58153(31) 

Freq 4th deriv, (10 -36 s -5 ) 

-0.690(14) 

Proj. semi-major axis, x (lt-s) 

8819(1250) 

11151(1250) 

Freq 5th deriv, i/g (10 -44 s -6 ) 

-1.26(9) 

Eccentricity, e 

0.93(3) 

0.95(3) 

Freq 6th deriv, z/g (10 -52 s -7 ) 
Freq 7th deriv, z/j, (10 -60 s -8 ) 

-6.2(3) 

-27(3) 

Long, of Periastron u (deg) 

12(4) 

10(4) 

Rms timing residual, a (ms) 

0.53 


0.57 

0.56 


data. A noise source was added to the data, defined by 
the power spectral density: 


P(f) =P™ + A 




( 3 ) 


where the white noise, P w = 6.5 x 10 24 yr 3 , and the 
power-law parameters are A = 9 x 10 _2 °yr 3 , f c = 
0.2 yr -1 and a = —4. The noise parameters were cho¬ 
sen to match the observed noise spectrum. The toasim 
plugins from tempo2 were used to generate 10 4 realisa¬ 
tions of the noise and the full fitting process was applied 
to each realisation. The error estimates were determined 
from the variance of the resultant fit parameters. 

Fig. le shows the timing residuals for the binary fit 
for / m = 10 Mq given in Table 1, column 4, showing 
that the rms is almost entirely limited by measurement 
errors, with any remaining unmodellcd systematic trends 
at a level that is comparable with those from pulsars of a 
similar age. We note that these binary models involve the 
fitting of 14 free parameters to achieve an rms residual 
of 0.57 ms. This is one parameter fewer than required for 
the timing-noise model with similar rms residuals that 
is presented in Fig. lc and Table 1, and based upon a 
7-derivative polynomial model, achieving an rms residual 
of 0.53 ms. Fig. Id shows the residuals based upon a 6- 
derivative timing-noise model that has the same number 
of fitted parameters as the binary model and is clearly a 
poorer description of the data. 

We conclude that there is a range of simple binary 


models that describe the TOA data very well. It is also 
clear that the reason we cannot define the orbit more 
precisely at present is that, because the data span covers 
only about 1/5 of an orbit, the differences between the 
models are buried within the measurement errors of the 
TOAs and likely timing noise. 

However, Figs. 3b-3f show that the values of many 
of the fitted parameters do not change much near the 
rms residual minimum. For all the models near this min¬ 
imum (Fig. 3a), the eccentricity increases from about 0.90 
to 0.95 with increasing / m (Fig. 3c). However, the corre¬ 
sponding To lies in a restricted range between MJD 58150 
and 58200 (Figs. 3b and 4). Thus the next periastron pas¬ 
sage will occur sometime during 2018 February or March. 
The intrinsic frequency derivative for all these models 
is about —6(1) X 10 -13 s~ 2 (Fig. 3d), corresponding to 
a characteristic age r = iz/2z> of about 180 kyr, which 
is substantially greater than the (rapidly decreasing!) 
characteristic age of the pulsar in the young-isolated- 
pulsar model, which has changed from ~ 120 kyr to 
~ 55 kyr over the six years. The values of right ascension 
a and declination S are not strongly model-dependent 
(Figs. 3e and 3f) and are about 20 h 32 m 13 s .106(8) and 
+41°27 , 24”.36(10), consistent with the position of the 
Be star, MT91 213. which is at 20 h 32 m 13 s .137( 18) and 
+41°27 , 24".48(20f]. an d which [Camilo et al.l (j2009l f con- 


1 Owing to a typographical error, ICamilo et al.l d2009h listed 
the last digits of declination of MT91 213 as 24.28; the correct 
value is 24.48 as used here. 
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Figure 3. The results of fits to the TOAs of a number of 
binary models having different fixed values of orbital period 
Pi j and mass function / m . a) The rms timing residual (cr) 
with the colour-coding the same as that defined in (c). b) The 
epoch of periastron (To), c) The orbital eccentricity (e). d) 
The pulsar rotational frequency derivative (i>). e) The right 
ascension of the pulsar less 20 h 32 m (A«). f) The declination 
less 41°27' (AS). 


sidered, but rejected, as a possible companion of the pul¬ 
sar. The conflicting conclusion of that study with this 
paper arose from their assumption of a circular orbit and 
the happenstance that their observations were conducted 
near apastron, where the gravitational influence of the 
companion on the pulsar is small, so that Doppler effects 
result in a small value of the pulsar second rotational 
frequency derivative. 
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Figure 4. The orbital variation of the frequency v of 
PSR J2032+4127 for the best fit binary models, for the in¬ 
dicated orbital periods Pu and having mass function / m = 10 
Mq. Note that the epochs of the next periastron around 58150 
are almost independent of the value of P b . The bold sections 
in the plots indicate the range of the available data. 

4 DISCUSSION 

If MT91 213 is indeed the companion star, what ad- 
dition al constraints might be placed upon the system? 
iKiminki et al.l (|2007h estimate the mass of the star to be 
M c ~ 17.5 Mg, based upon their estimate of its spec - 
tral type (BOV) and the study by iMartins et al.l (120051) 
who suggest that such spectrally-determined mass esti- 
mates may b e in er ror by 35-50%. However, more recently 
iHohle et al.l (|2010l ) have estimated the mass of BOV stars 
to be 15.0±2.8 Mq using 2MASS spectrometry and Hip- 
parchos parallax data. We adopt this latter value in the 
following discussion. Unfortunately, the inclination of the 
plane of the orbit to the plane of the sky is unknown. 
For random orientations of a binary orbit, the median of 
the probability density function of inclination i occurs at 
% = 60°. For an orbit with this inclination and a pulsar 
of mass M p = 1.35 Mq, the mass function (Equation [lj 
would thus be / m = 8.2 Mq, which corresponds approx¬ 
imately to the first binary model in Table 1, column 4, 
and the orange lines in Fig. 3. The greatest possible value 
of the mass function for these two stars is about 15 Mq 
which occurs if i = 90°. There is no doubt that an orbit 
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with a mass function of ~ 5 — 15 Mg, having an orbital 
period lying between 7000 and 11000 days, or about 20- 
30 years, is entirely consistent with the data (Fig. 4). 
Although larger mass functions also provide satisfactory 
fits to the data, they are inconsistent with a companion 
that has the mass of MT91 213. 

Two other pulsar/Be-st ar binary system s are k nown, 
PSR B1259—63/S S 2883 lljoh nston et al.1 Il992h and 
PSR J1638—4725 llLorimer et al. 20061 ; Lvnal200Sl) . hav¬ 
ing orbital periods Pb of 1237 d and 1941 d respectively, 
and companion masses M c of ~ 4 Mg and ~ 8 Mg. 
The PSR J2032+4127/MT91 213 system is more extreme 
than either of these, with Pb ~ 9000 d and M c ~ 15 Mg. 
Fig. 5 shows the most likely orbital configuration. The 
pulsar is presently on the Earth side of the Be star and 
is beginning to move rapidly in towards periastron, after 
which it will move behind the star and any circumstel- 
lar disk. Such stars have substantial stellar winds and, in 
the cases of both PSR B1259—63 and PSR J1638—4725, 
as well as radio eclipses, there are increases in both dis¬ 
persion measure (DM) and multi-path scattering close to 
periastron, which arise as the pulsar becomes more em¬ 
bedded in the dense circumstellar environment. No sig¬ 
nificant change has yet been detected in the scattering 
or DM of PSR J2032+4127, with the fitted values of its 
first derivative DM1 consistent with zero (Table 1), but 
it is still far from periastron and this is not unexpected. 
Following periastron in early 2018, the pulsar will move 
behind MT91 213 and may be eclipsed, at least in the 
radio, by its atmosphere or disk, depending upon the in¬ 
clination of the orbit to the line-of-sight. 

Such a large, long-period system is unique, but its ex¬ 
istence and its survival of the velocity kick that the pulsar 
probably experienced when it was formed is perhaps not 
surprising, because of the deep gravitational well of the 
~ 15-Mq companion. The separation of the markers in 
Fig. 5 gives an indication of the orbital velocity, showing 
that, although the velocity at apastron is only a few km 
s -1 , at periastron the velocity is in excess of 100 km s^ 1 . 
It is a system that, unlike the majority of pulsars, did not 
quite become unbound after the neutron star formation. 

In all the allowed binary models, periastron will oc¬ 
cur in just over three years time. Much of the degeneracy 
present in the fits will be resolved as periastron is ap¬ 
proached. In particular, the rapid increase in slow-down 
rate during the next 1-2 years will establish the mag¬ 
nitude of the eccentricity and hence the stellar separa¬ 
tion at periastron. All evidence indicates that MT91 213 
is in the Cyg OB2 stellar association, which is located 
at a distance of 1.4-1.7 kpc, determined both by spec¬ 
troscopic parallax and by VLBI determination of the 
trigonometric parallax of masers associated with t he as¬ 
sociation stars (fjVIassev fe Thompsonlll99ll : lHanso nll2003l: 
iRygl et a,l.ll2012t TT in the binary hypothesis, this is there¬ 

fore also the distance of PSR J2032+4127. We note that, 
at this distance, the projected major axis of the pulsar 
binary orbit will subtend about 25 mas viewed from the 
Earth. It should be possible to track the pulsar around 
the orbit using a VLBI network such as the VLBA or 
the EVN, hence establishing the size and the projected 
eccentricity of the orbit. Together with timing measure¬ 
ments, such observations permit the direct determination 


of the inclination i of the orbit to the plane of the sky 
and the distance of the system from the Earth. 

The spin-down luminosity of PSR J2032+4127 in the 
binary model s presented here is E = 1.7 x 10 35 ergs _1 , 
62% of the lCamilo et al] (l2009ll isolated-pulsar value. The 
efficiency in converting rotational kinetic energy to 7 -ray 
luminosity P 7 , 7 = L^/E, is therefore increased by a fac¬ 
tor of 1.6. Using the distan ce of 1.4 1.7 kpc and the up¬ 
dated L 7 from lAbdo et al.l (120131 4. and assuming as usual 
an isotropic 7 -ray beam, then 7 =15%-20%, which is un¬ 
remarkable for a young puls a r (t he much higher value 
of 77 presented in lAbdo et al.l (j2013l ~) uses the DM-based 
distance of 3.7kpc). 

We also note that this system is still fairly young, 
the pulsar having a characteristic age of 180 kyr and the 
Be star having a total lifetime of only a few million years, 
which needs to encompass the lifetimes of the pulsar pro¬ 
genitor star and of the pulsar itself. 

The fact that PSR J2032+4127 is a member of a 
long period binary system with a massive companion star 
means that we should revisit the interpretation of the spa- 
tially coincident e xtended TeV 7 -rav and X-ray emission 
(I Aliia et al.ll2014h . iKargaltsev et alT (120141 ') have recently 
shown that the similar system PSR B1259—63 exhibits 
variable extended X-ray emission that is located outside 
the binary system and is moving away. It may be that this 
is emission that is associated with a circum-binary shock 
which can be gener ated between the wi nds of a pulsar 
and a massive star dBosch-Ramon et al.ll20l3) . We note 
also that if the interpretation of the TeV emission exten¬ 
sion is to do with the proper motion of the binary then 
the increased age will result in a reduced transverse ve- 
locity_of less than 30 kms -1 compared to that derived by 
lAliu et ail (l2014f) . 


In conclusion, there is a restricted range of binary or¬ 
bits that describe the observed timing data very well and 
explain the extraordinary and unique growth in the ob¬ 
served rate of rotational slow-down. This and the coinci¬ 
dence of the precise positions of the pulsar and MT91 213 
all point to the two stars comprising a neutron-star/Be- 
star binary system in the Cyg OB2 association. Near peri¬ 
astron, the pulsar may well become obscured in the radio 
by free-free absorption and by severe pulse scattering in 
the Be-star circumstellar wind and disk as it passes be¬ 
hind the star, in the same way as the Be star SS 2883 
obscu res PSR B1259—63 at periastron lljohnston et al.l 
1 19961 4 . However, studies of PSR J2032+4127’s DM and 
rotation measure (RM) variations are likely to present a 
rare opportunity to study the density of the stellar wind 
and any circumstellar disk and the magnetic field of the 
Be star. Although it may become obscured in the radio 
for a short while, it should be possible to track its rota¬ 
tion and orbit around periastron by observation in 7 rays. 
Close to periastron, the MT91 213 Be star may also dis¬ 
play optical and/or X-ray variability as a direct result of 
its interaction with the pul sar. Even now MT91 213 dis¬ 
plays X-ray vari ability fe.g.lR auw et al.ll2014ll and optical 
varia bility (e.g. ISalas et all 20141 ; see also Camilo et al.l 
120091 1 . although these are not atypical of the intrinsic 
variability observed in other Be stars and almost surely 
have nothing to do with the neutron star. 
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Figure 5. Schematic diagram showing the approximate orbital motions of PSR J2032+4127 and its Be-star companion MT91 213 
about their common centre of mass for the binary model given in Table 1, column 4 for / m = 10 Mq, projected onto the plane 
containing the line-of-sight and the major axis of the orbit. The inclination i of the plane of the orbit to the plane of the sky is 
assumed to be 60°. The markers are at 200-day intervals and indicate the time from periastron. The thick line shows the portion of 
the pulsar orbit covered by the observations reported here. The pulsar is moving counter-clockwise in the diagram and is currently 
(MJD 56912) on the near side of the Be star and about 1300 days before periastron passage. Note that the orbital velocity is 
proportional to the separation of the markers, a 1000-light-sec separation indicating a velocity of about 18 km s —1 . The small 
ellipse near the origin shows the orbit of the Be star, assuming that it has a mass of 15 Mq and that the pulsar has a mass of 
1.35 Mq. At an estimated distance of 1.7 kpc, the projected major axis of the pulsar orbit will subtend an angle of about 25 mas. 
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